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We examine the effect of d-wave symmetry on zero bias anomalies in normal-superconducting 
tunnel junctions and phase-periodic conductances in Andreev interferometers. In the presence of 
d-wave pairing, zero-bias anomalies are suppressed compared with the s-wave case. For Andreev 
interferometers with aligned islands, the phase-periodic conductance is insensistive to the nature of 
the pairing, whereas for non-aligned islands, the nature of zero-phase extremum is reversed. 
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■ During the past few years, studies of the subgap conductance of normal-insulator-superconductor (N-I-S) structures 
have revealed new and unexpected behaviour Whereas conventional tunneling theory suggests that the subgap 

■ conductance must vanish, experiments reveal that when the normal region becomes phase-coherent, there exists a 
zero- voltage peak in the differential conductance, which can be comparable with the normal state value. The interplay 
between Andreev scattering j2| at an N-S interface and disorder-induced scattering in the normal region has been 

| recognised as the main physical origin of this zero bias anomaly (ZBA) ||- pjj ]. 

Andreev scattering |2| involves the simultaneous tunneling of two electrons of opposite momenta through the N-I-S 
structure. In the case of a low-conductance tunnel junction (G tun <C 1), this two-particle process occurs with a 
probability « G\ un . As a consequence, the conductance in the superconducting state is expected to be much smaller 
than that in the normal state. However, disorder-induced scattering on the normal side may give rise to particle- 
particle correlations, which effectively, increase the probability for two particles with opposite momenta to tunnel 
into the superconductor. This increase manifests itself as a low total- momentum singularity in the particle-particle 
I i scattering channel (the so-called Cooperon) jljfl , which is also the relevant scattering channel for the occurrence of 
s-wave pairing in BCS theory jlqj . In the latter case, the low total- momentum instability signals the formation of 
a bound state of two electrons with opposite momenta, while in the former it leads to enhanced backscattering in 
disordered systems. In both cases, the singularity in the particle-particle Cooperon channel is characterized by the 
s-wave symmetry of the relative two-particle wavefunction and is ultimately responsible for the occurrence of the ZBA 
||. In contrast, for a clean normal region, Andreev scattering is almost suppressed |14| and conventional tunneling 
theory applies. 

The aim of this paper is to demonstrate that ZBAs and related phenomena are sensitive to the symmetry of the 
superconducting order parameter and therefore can be used to distinguish between s-wave and d-wave pairing in high- 
T c superconductors. In what follows, we extend a numerical multiple scattering approach |l5[ | to the case of non-local 
anisotropic superconductors. In the case of local s-wave pairing, this approach has been shown to be equivalent to 
other techniques and has been used to study the cross-over between different transport regimes Jl6[] . To start with, 
in section II , we present an analytic treatment of a ballistic N-I-S structure, which generalizes the theory of Ref. 
Jblj to the case of a tight-binding lattice and of a non-local superconducting order parameter. In section III, we 
show that, in the presence of disorder, the subgap conductance of a normal-insulating- d-wave junction is suppressed 
compared with the corresponding normal-insulating-s-wave structure. In particular, in the regime of small tunnel 
junction conductance Gtun <C 1, we predict that the d-wave junction shows only a weak zero bias anomaly. Similarly, 
the zero energy dip of a N-I-S structure, in the high tunnel junction conductance regime Gtun ~ 1, is considerably 
smaller for the d-wave case. 

Having examined the case of a single superconducting contact, in section IV, we extend the analysis to a phase- 
coherent structure in contact with two superconductors. The electrical conductance of such structures is known to 
be a periodic function of the difference between the order parameters phases of the two superconductors, which in 
turn is an externally controllable quantity. Such Andreev interferometers have recently been the subject of intensive 
theoretical ]r7|| - [24J and experimental studies J25|| - 28 1. In what follows, we predict that certain features of the 
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phase-periodic conductance, such as the nature of the zero-phase extremum, are sensistive to the symmetry of the 
order parameter, and, therefore, interferometers of this kind provide a further probe into the nature of the pairing. 



II. ANDREEV REFLECTION AT A N-S INTERFACE 



To begin with, in this section, we examine a normal-superconducting interface (N-S) in the presence of a non-local 
pairing potential, described by the Hamiltonian [ p9[ 

H = Yl £ * C L C ^ + Y [ic{<y c i+&,<r + h - C ] + Y [( C h4+S,l ~ C l,l4+S,f) A ^,S + h ' C - ' W 

Here, the index i runs over a two-dimensional tight-binding lattice, with unit lattice constant, 8 sums over nearest 
neighbours, 5 = x,y, and the pairing potential, Aj )( 5, is defined on the bond from site i to site i + 5. The operator 
cj CT {ci,a ) creates (destroys) an electron at site i, with site energy ej, and 7 is the hopping matrix element between 
nearest-neighbours sites. All energies will be measured in units of 7, which will be set to unity throughout the paper. 
The above Hamiltonian is diagonalised by solving the Bogoliubov - de Gennes equation: 

Efa = -Cifa - Y,g 7 {<t>i+s + <Pi-s) + ( A i,W>i+<5 + Af-^i-i) i U1 

where tpi(4>i) indicates the particle (hole) wavefunction. In general, the Hamiltonian of eq.(Q) is a mean-field ap- 
proximation to a more complex Hamiltonian containing electron-electron interactions and all parameters should be 
determined self-consistently. However, in many cases of experimental interest, the qualitative form of parameters such 
as Aj j( 5 is known and for the purpose of highlighting generic transport properties, self-consistency is not required [j30| . 

The Bogoliubov-de Gennes equation may be solved by means of a transfer matrix method or a recursive Green 
function technique. These methods work in any dimension and constitute the only exact approach when translational 
invariance is absent. As a prelude to such a calculation, in this section, we begin by considering a system with 
translational invariance in the direction perpendicular to the current flow. Our motivation for doing so is two-fold. 
On the one hand, one obtains new results for Andreev scattering in the presence of anisotropic pairing. On the other, 
the analysis provides a controllable limit, which can be used to test the numerical machinery used in more complicated 
situations. 

In the presence of translational invariance in the direction transverse to the current, the problem reduces to one 
of many independent one-dimensional channels, each characterised by one or more discrete quantum numbers. In 
what follows, we consider the case of a 2-dimensional system with a N-S interface, whose normal vector points in 
the x direction, although by redefining the parameters p,, A introduced below, all results are trivially generalized to 
3-dimensions. The number of independent channels is determined by the width M of the system, and each channel 
has a discrete wavevector k y , along the y direction. Choosing ej = — fj, yields an energy dispersion relation in the 
normal region of the form 

e k = -2"/(cos(k x ) + cos{k y )) - fi (3) 

where the uniform site energy /i determines the filling of the tight-binding band. Writing 

p, = [i + 2"fcos{k y ). (4) 

shows that the dispersion relation reduces to that of a one-dimensional system with a channel-dependent chemical 
potential p. 

In the superconducting region, the dispersion relation (||) is replaced by 

E 2 = (e- k f + \A- k \ 2 , (5) 
where e k — —2 r y(cos(k x ) + cos(k y )) — fi, and the momentum dependent gap function, Ajr,, is given by 

Aj. = 2(A x cos(k x ) + A y cos(ky)), (6) 

where for d-wave symmetry, A x — —A y . 
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In writing the above equation, we have assumed a uniform pairing potential for all sites i > 0. The transverse 
component, k y is conserved through the interface, k y = k y , and for a fixed value of k y , the gap function can be written 

A- k = 2A x cos{k x ) + A, (7) 

where A = 2A y cos(k y ), which demonstrates that the pairing potential along the y direction yields a local contribution 
to the pairing potential for a particular channel. Hence for i < 0, the Bogoliubov - de Gennes equation reduces to 



Eipi = -7(^-1 + V'i+i) ~ WAi 
E(j)i = 7(^-1 + <fii+i) +Ji4>i 



(8) 



while for i > 0, 



Eipi = -7(^-1 + V'i+i) - fJ'i'i + ^4>i + &x(<f>i-l + <f>i+l) /qx 

E<f>i = 7(^_i+^ i+ i)+^ + A^ j + A !B (^ i _i+Vi+i). 

Since the value of the transverse component of the momentum enters only through the effective chemical potential 
p and the local contribution to the gap A, it is convenient to drop the x suffix in labelling the various longitudinal 
momenta and write k and q, (k and q), for the longitudinal momenta of particles and holes in the normal (super- 
conducting) region. With this notation and writing (ip R ,<f) R ) and (ip L ,(f) L ) for the solutions when i > and i < 0, 
respectively, one obtains for i < 

I L _ ikRi 1 -ikRi 

and for i > 0, 

The coherence factors u and t; identify a particle-like excitation of energy E and momentum fc, while u and w 
correspond to an hole-like excitation at the same energy and momentum q. To compute these quantities one notes 
that 

E 2 = (-2 7 cos(p) - p) 2 + |A + 2A x cos(p)\ 2 (12) 
where p may be k or q, which yields 

cos{p) = -\{p + A±B)/ 1 (13) 

with 



A = A. 



7 A - A^ 

: 7 2 + A^ ' 



S = 7,, 



£ 2 f 7A-A a /2 
^7 2 + A 2 y 7 2 + A 2 ^ 



For the coherence factors one has to distinguish the case of real B, which corresponds to quasi-particle transmission 
through the interface, from the case of imaginary B, which corresponds to no quasi-particle transmission. 
For B real, one obtains 



u 2 (v 2 ) 
u 2 {v 2 ) 



Ul + {-){A + B)/E) 
\(l + (-){A-B)/E) 



(14) 



with u/v = sign {Aj./{E — ej)), and u/v = sign (Ag/(E — 6g)). 
For B imaginary, we write B = iB, to yield 

u/v = Aj t /(E -A- iB) 
u/v = Ag/(E -A + iB) = (u/v)* 



(15) 
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To solve for the scattering coefficients r , r a , t , and t a , one needs matching conditions at the interface. These are 
obtained by evaluating equations (^) and (Q) at i = and at i = — 1 and can be written in the form 

= <p« (16) 

and 

<i> L - 1 =<t>% (17) 

where (see appendix) 4> I ^ 1 , are obtained by acting on <f>f-, ip^ with the appropriate transfer matrix at the interface. 
Equation ( |l6| ) yields 

1 + r a = t a u + t a u . . 

r a = t D v + t a V, 

whereas the matching conditions eq.([r^) yield (see appendix) 

e- tk + r a e ik = t ue~ ik (l - (w/«)(A B / 7 )) + t a ue^(l - (v/u)(A x /j)) 

r a e-*« = t ve- lk (l + (u/v)(A x /~/)) + t a ve^(l + (u/v)(A x / 7 )). [ ' 

The eqs.( ^8| - |l9| ) yield all the scattering coefficients associated with a particle incident from the left on a clean N-S 
interface. A direct analytic evaluation is rather messy, but it is trivial to solve these numerically. Explicit results will 
be presented in section III. In the presence of a tunnel barrier modelled by a delta function potential at the interface, 
ti is replaced by ej = — /z + U5i t o and the matching conditions eq.(|l9|) by 



e - lfc + roCik = tolte -ik(i _ ( w / u )(A x / 7 )) + t a ue^(l - (v/u)(A x / 7 )) + (1 + v )Uh 
r a e-^ = t a ve- i% (l + (w/u)(A x / 7 )) + t a ve^(l + (m/u)(A x / 7 )) + r a U/j. 



Various limiting forms of the above expressions are discussed in appendix. Here we merely note that with the 
convention adopted in eqs.( ^0[]Tl| ), if Vk (v q ) and vj, (vq) are group velocities for particles (holes) in the normal and 
superconducting regions respectively, then the following unitarity condition is satisfied, R a + Rq + T a + To = 1, where 
Ra = (v a /v k )\r a \ 2 , R a = |r | 2 , T a = {v q /v k )\t a \ 2 and T = {v k /v k )\U\ 2 . 



III. RESULTS FOR N-I-S STRUCTURES 

In this section we present explicit results for the above scattering coefficients and for the electrical conductance 

G = (x) ^ (1 ~ Ro + Ra) (21) 

ky 

where the sum is over all channels. 

For an N-S interface with no potential barrier, fig.|l|(a) shows the behaviour of the differential conductance as 
a function of energy, along with various scattering coefficients. In a conventional superconductor under sub-gap 
conditions, only Andreev reflection contributes to the current flow, because sub-gap quasi-particle transmission is 
forbidden. In a gapless superconductor, the situation is more complicated, because each channel has its own effective 
gap, so that normal quasi-particle transmission occurs even at low energies. To illustrate this point, fig.^(a) shows the 
behaviour of transmission and reflection probabilities for a system width M ~ 10. The number of open channels N 
depends on the position of the chemical potential fj, within the band. Here we have used fi = — O.27 so that N — 9. For 
free-end boundary conditions in the transverse direction, the allowed values of k y are nn/ (M + I), with n = 1, 2...M. 
For this choice, all channels have non-vanishing gap and at zero energy, normal transmission vanishes. By increasing 
the energy, one eventually crosses the effective gap of a particular channel, at which point transmitting channels 
appear in the superconducting region and Andreev scattering is suppressed. Since Andreev scattering contributes a 
factor of 2 (in units of 2e 2 /h) to the electrical conductance, while transmission processes only contribute a factor of 
unity, the conductance decreases at such energies. In fig.|l|(b), we show results for a system width M = 1000. At zero 
energy, it is interesting to note that the conductance per channel is insensitive to the width of the system and agrees 
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perfectly with the M — 10 value. At finite energy the results differ slightly, because in the case of a large number of 
channels, channels open continously with increasing energy, to yield the smooth behaviour shown in fig.[j](b). 

We now consider the case of an N-I-S structure, in which a tunnel junction is present at the interface. Fig.||(a-b) 
shows results for the electrical conductance G and for various scattering coefficients in the presence of a barrier height 
U = 2.O7 for a system width M = 10. The oscillating behavior is a finite size effect which can be understood by 
observing that a peak occurs when the energy becomes greater than the effective gap A of a particular channel. The 
oscillation arises because for E < A the contribution from such a channel is proportional to the square of the barrier 
transmission coefficient T (i.e., w G 2 un ), whereas for E > A it is proportional to the first power of T (i.e., » G t „„). 

In fig.|^(c-d) we show the conductance for a width M = 1000. Again the conductance per channel, at zero energy, is 
found to be insensitive to the number of open channels, and oscillations at intermediate energies are no-longer present. 

Having examined a clean metal in contact with a superconductor, we now introduce disorder to the normal metal, 
to produce a diffusive conductor in contact with a d-wave superconductor. In this case the conductance is obtained 
using the transfer matrix method outlined in reference fl32"f and computer resources impose restrictions on the system 
width M. 

The N-I-S structure of fig.§ consists of a diffusive metallic region placed in series with a tunnel junction, which 
in turn is adjacent to a superconductor. In the absence of disorder, the numerical code agrees exactly with the 
analytical results of figures |l| and ||. In what follows, the simulated structure is a two-dimensional tight-binding 
lattice of width M = 10 sites. The disordered region is of length Ldif sites, the tunnel junction is L tun sites long 
and the superconductor has a length L sup . The conductance of the entire structure is denoted by G and the average 
over an ensemble of disorder realizations is < G >. The physical variables in the following calculation are the 
averaged conductance < Gdif > of the diffusive region and the conductance of the tunnel junction Gtun- To identify a 
suitable choice of parameters, we considered first a normal diffusive portion of length Ldif and width M , connected to 
crystalline, normal leads. The conductance of a diffusive material is inversely proportional to its length and therefore 
a plot of < Gdif > Ldif as a function of Ldif will exhibit a plateau in the diffusive regime, with a mean free path 
given by I =< Gdif > Ldif/M(2e 2 /h). A diffusive system must satisfy I <C Ldif an d I *C M. Furthermore, if weak 
localization corrections are to be neglected, we require Nl 3> Ldif. In the calculations which follow, having in mind 
also the necessity of minimizing the CPU time, we have made the following choice of parameters: Ldif — 30 sites, 
disorder width W — 1, system width M = 10, \i = so that the number of open channels is N = 10. This yields a 
mean free path I = 4.2 and an average conductance < Gdif >= 1-6 (in units of 2e 2 /h). The superconductor has a 
length L sup = 100 with order parameter A x = O.I7. One characteristic energy scale is the maximum gap, which for 
A x = O.I7 is equal to OT7. The tunnel junction is L tun — 1 site long and the potential on the line of sites defining 
the junction is — ^ + with eb taking the value 77 and 27, which yields Gtun = 0.4(2e 2 //i) and Gtun = 3(2e 2 //i) in 
the low and high tunnel junction regimes, respectively. 

For the case of low tunnel junction conductance, fig.(^) shows the sub-gap conductance as a function of the energy, 
which in the linear response regime, corresponds to experimentally measured I-V characteristic. The strong peak, 
present in the s-wave case, is strongly suppressed in the d-wave case. Fig.(||) shows the corresponding behaviour in 
the high tunnel junction conductance regime. In this case, the zero energy behaviour is characterized by a dip in the 
conductance. As before, the d-wave conductance and the zero bias feature is suppressed. 



IV. RESULTS FOR D-WAVE ANDREEV INTERFEROMETERS. 



We now consider the effect of d-wave symmetry on the properties of Andreev interferometers. Two different 
interferometer geometries are analyzed below and shown as inserts in figures 6 and 7. In each example, the system 
consists of two superconducting regions with order parameter phases <pi and <p2, separated by a normal region TV, 
with a quasi-particle current flowing vertically. In figure 6, the S-N-S structure is placed in contact with a tunnel 
junction and the N and S regions are clean. In figure 7, there is no tunnel junction, but the whole S-N-S structure is 
disordered. 

The physical parameters used to obtain figure 6 are as follows: fi = and the tunnel barrier site-energy e& = 27, 
which corresponds to an average conductance per channel of Gtun/N = 0.176(2e 2 //i). As shown in figure 6, the 
electrical conductance in the presence of s-wave pairing shows a large amplitude of oscillation and a zero-phase 
minimum, which are characteristic of a ballistic structure, as discussed in Refs. Jl6| , p3j . In contrast, for a d-wave 
interferometer, the nature of the zero phase extremum depends on the relative orientation of the two superconducting 
islands. In one case (solid line in figure ^) the x axes with a positive value of the order parameter are parallel in both 
islands, while in the other case (dashed line in figure ^|) the x and y axes are exchanged in going from one island to the 
other. For the case of parallel orientation, the d-wave interferometer is almost identical to the s-wave case, exhibiting 
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a large amplitude of oscillation and a zero-phase minimum. When the islands are orientated perpendicular to one 
another, the minimum becomes a maximum and the entire curve is shifted by n. This arises because an electron 
Andreev reflected at an N-S interface, aquires the phase of the bonds in the longitudinal direction (with respect to 
the incoming direction of the electron). When the two islands are orientated parallel, they have the same value for 
the longitudinal bonds, so that there is no effective phase difference between the two islands, as in the s-wave case. 
When one of the two superconducting islands is rotated by ir/2, the longitudinal bonds have a phase difference of w, 
which adds to the external phase difference. This effect is similar in origin to that predicted |53| and later verified |34|] 
in a corner SQUID experiment on the high-T c superconducting compound YBCO. The relevance of the orientation of 
the crystal axis with respect to the N-S interface has been also recently put forward |35|] as a possible explanation of 
zero bias anomalies. 

For figure 7 we choose parameters corresponding to a diffusive region. The structure has a length L = 50 with 
each of the three regions having a width M' = 15, so that the entire structure has a width M = 45. We used /i = 0, 
disorder width W = O.57, which yields, in the normal state, a mean free path I = 12. For an s-wave interferometer 
(shown as an inset), the zero phase extremum switches to a maximum and the value at <j) = ~n becomes a minimum. 
For aligned islands (solid line), the d-wave interferometer shows the same qualitative behaviour as in the the s-wave 
case. However, rotating one of the islands by 7r/2 again shifts the curve by 7r and changes the zero-phase extremum 
from a maximum to a minimum. 



V. DISCUSSION 



Zero bias anomalies and phase-periodic conductances in Andreev interferometers are paradigms of phase-coherent 
transport in hybrid N-S structures. In this paper we have examined the sensistivity of these phenomena to the 
symmetry of the superconducting order parameter. We find that ZBAs are suppressed, reflecting the fact that 
disorder induced scattering and Cooper pairing no longer occur in the same particle-particle scattering channel. We 
also find that for aligned islands, Andreev interferometers are relatively insensistive to the nature of the pairing. 
Neverthless, such devices could be used to reveal the presence of d-wave pairing, since the positions of conductance 
extrema are sensitive to the relative orientation of the two order parameters. 

For simplicity, we have avoided the necessity of a fully self-consistent calculation, by restricting the analysis to a 
particular crystal orientation. Analyses |30| of the proximity effect in ballistic systems reveal that a self-consistent 
theory can yield a drastic suppression of the superconducting order parameter in the vicinity of the N-I-S interface, 
but only for certain crystal orientations of the d-wave superconducting order parameter with respect to the interface 
p6| . In this paper, our main focus has been the interplay of Andreev scattering at the interface and disorder- induced 
scattering in the normal region, and therefore we have not considered these particular crystal orientations. For this 
reason we expect that our main conclusions will not be qualitatively changed by a fully self-consistent calculation, 
but for the future it would be of interest to explore transport in the presence of alternative crystal orientations. 
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APPENDIX A: TECHNICAL DETAILS OF THE ANALYSIS OF SECTION II 

To obtain equation (|l^) one notes that the left-hand sides of eq.(|l^) are 

fa = -((E + p)/i)fa - fa = e- lk + r a e lk 
fa = {{E - Ji)h)fa - <h L ~ = r-«-<« 



whereas the right-hand sides are 
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ifa = -(E + p,)il>g - 70f + &</>§■ + A a ^f 
7 fl ! ={E- fittf - 7 < - Aip§ - A^f . 

Notice that i/i^ and 0^ are decoupled because there is no superconducting bond between i = — 1 and i = 0. 
Inserting eq.(|ll]) into the above expressions and using the eigenvalue equations, yields 

fa = t ue-^{l - {v/u){A x /j)) +t a ue^{\ - (v/u){A x / 7 )) 
fa = t ve- lk {\ + (u/v)(A x / 7 )) +t a ve^(l + (u/v)(A x /~/)) 



Finally one obtains the last two matching conditions in the form (|l9|). 

The equations obtained in section II are very general and yield a variety of useful results. To illustrate this we first 
rewrite equation (pfl) in the form 



f e' fc - e -"°(l-(t-/ M )(A x / 7 ))-(j7/ 7 ) , e» fc - e "?(l-(e/fi)(A x / 7 )MC// 7 ) _ , 
l oU e ik_ e -ik r T, a U e ik_ e -ik — J- 

e -^-e-^(l+( M / t ,)(A x / 7 ))-(i7/ 7 ) e-^-e^(l + (fi/e)(A 3; /7))-(t//7) _ n 

l O» „sfc _„-ifc T to" „ik_„-ik — <J. 



(Al) 



Once the above pair of equations is solved for t Q and t a , eq.(|Tq) yields r and r a . In the case of no quasi-particle 



transmission, the only physical quantities are r Q and r a , which, from the form of eq.(Al), depend only upon the ratios 
u/v and u/v. 

Results corresponding to the case of a local order parameter A Q can be obtained from the above eqs.(]Al|) by simply 
setting A x = 0, and A = A G . By considering for simplicity an ideal interface (U — 0), one obtains 



t Q u[ 



t V e Jk_^-ik +t a V\ 



(A2) 



which gives 



t = v(e^ 



k )/d 



(A3) 



and, by using eq.(]Ig|), 



t a = v{e- iq - e~ lk )(e lk - e~ tk )/d 



r Q = \vu(e is - e- iq ){e- fk - e - lk ) + uv( e - lq - e- fk )(e lq - e^)! /d 



(A4) 



(A5) 



vv(e 



ig 



-ik\ ( Ah 



- -i k 



)/d 



(A6) 



where 



d = 



vu(e zk - e~ zk ){e lq - e^) - uv(e lk - e zq )( 



e- lk - e 



Eqs.(|A3|- A6) solve the problem of determining the various scattering coefficients. As a last step one has to substitute 
the expressions for the momenta k, q, k, and q in terms of the energy E as given by eq. (|l3|) . We will not do this 
substitution here, because the resulting equations look rather cumbersome. We will instead content ourselves by 
considering the relevant limit of zero quasiparticle energy. In this case there is no quasiparticle transmission through 
the interface, and the momenta k and q are complex. By setting k{q) = p + (— )il and observing that q — k and 
j2 = —2cos(k) = —2cos(p)cosh(l), Aq = 2sin(p)sinh(l) yields for the Andreev reflection scattering probability 



R a = {v q /v q )\r a \ 2 = 



2(1 - OiW)) 



1 - (/2 2 /47 2 ) + (A2/4 7 2) + + (/22 + A2)/4 7 2))2 - £2 



(A7) 
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It is clear from the above equation that in the limit A Q — > 0, the Andreev scattering probability R a — > 1, which 
amounts to the so-called Andreev's approximation largely used in a number of theoretical treatments. To close this 



appendix, it is useful to make connection with the continuum limit, which is most conveniently obtained from eqs.(A3 



A6). To make this connection, we restore the lattice step a in all expressions involving the momenta, i.e., k — > ka and 
take the limit a — ► 0, 7 — > 00, such that 7a 2 = h 2 /2. This leads to 

r a = [uv(q + q){k-k)- uv(q -k)(q + k)] jd (A8) 

r a = [2vvk(k + q)] jd (A9) 

t a = [2vk(q + q)] jd (A10) 

t a = [-2vk(q-k)] /d (All) 

where 

d = uv{k + k)(q + q) — uv(q — k)(k — q). 

In the limit (A G //x -C 1), Andreev's approximation of ignoring both Andreev transmission and normal reflection is 
valid and one obtains [HJ] 

i? a = 1, E < A (A12) 

and 



^WF=1' s>Ao (A13) 

In contrast, at zero energy, without invoking Andreev's approximation, one obtains the exact result 

R a = ^_ (A14) 

with i? = 1 — i? a . The above equation clearly shows how Andreev's approximation breaks-down in the large A D 
regime. 
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FIG. 1. Electrical conductance (thick solid line), Andreev reflection (solid line line), normal transmission (dashed line). 
Normal reflection and Andreev transmission are negligible on this scale. Parameters used are: (a) width M — 10, number of 
open channels N = 9; (b) width M — 1000, number of open channels N = 857; in all cases fi — — O.27, A x = O.OI7, barrier 
height (7 = 0. 

FIG. 2. Electrical conductance (thick solid line), Andreev reflection (solid line), normal transmission (dashed line), figures 

(a) and (c); normal reflection (solid line), Andreev transmission (dashed line), figures (b) and (d). Parameters used: (a) and 

(b) width M = 10, number of open channels N = 9; (c) and (d) width M = 1000, number of open channels N = 857; in 
all cases n = — O.27, A^ = O.OI7, barrier height U — 27. Values are normalised to the normal state conductance in order to 
compare with the case of no potential barrier at the interface. 

FIG. 3. Schematic picture of an N-I-S stucture, considered in the text. Parameters used are: width M = 10, lenght of the 
diffusive region Ldiff = 30, length of the tunnel junction Lt un = 1, disorder width W = 1, fx = 0, number of open channels 
N = 10. 

FIG. 4. Total conductance as function of energy in the low tunnel junction conductance regime, s-wave: dotted line; d-wave: 
solid line. Gtun = 0.4 in units of 2e 2 /h. The conductance is normalised to the normal state conductance G . The two insets (left 
d-wave and right s-wave) show the conductance behaviour on a more extended energy range. The value of the superconducting 
order parameter in the s-wave case is A s = O.I7. 

FIG. 5. Total conductance as function of the energy in high tunnel junction conductance regime. S-wave: dotted line; d-wave: 
solid line. Gtun — 3 in units of 2e 2 fh. The conductance is normalised to the normal state conductance G . The two insets (left 
d-wave and right s-wave) show the conductance behaviour on a more extended energy range. The value of the superconducting 
order parameter in the s-wave case is A s = O.I7. 

FIG. 6. The figure shows results for a system with an insulating barrier placed at one end of the interferometer (See the 
structure displayed in the top right corner). The solid line (dashed line) indicates results for the conductance with aligned 
(not aligned) superconducting islands, for a d-wave interferometer. The inset shows the corresponding result for an s-wave 
interferometer. 
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FIG. 7. The figure shows the conductance for a system without a barrier but when the interferomter is dirty (See the structure 
displayed in the top right corner) . The solid line (dashed line) indicates results for the conductance with aligned (not aligned) 
superconducting islands, for a d-wave interferometer. The inset shows the corresponding result for an s-wave interferometer. 
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